RESURGENCE OF THE EULER-MACLAURIN SUMMATION FORMULA 



OVIDIU COSTIN AND STAVROS GAROUFALIDIS 



Abstract. The Euler-MacLaurin summation formula relates a sum of a function to a corresponding integral, 
with a remainder term. The remainder term has an asymptotic expansion, and for a typical analytic function, 
it is a divergent (Gevrey-1) series. Under some decay assumptions of the function in a half-plane (resp. in 
the vertical strip containing the summation interval), Hardy (resp. Abel-Plana) prove that the asymptotic 
expansion is a Borel summable series, and give an exact Euler-MacLaurin summation formula. 

Using a mild resurgence hypothesis for the function to be summed, we give a Borel summable transseries 
expression for the remainder term, as well as a Laplace integral formula, with an explicit integrand which 
is a resurgent function itself. In particular, our summation formula allows for resurgent functions with 
singularities in the vertical strip containing the summation interval. 

Finally, we give two applications of our results. One concerns the construction of solutions of linear 
difference equations with a small parameter. And another concerns resurgence of 1-dimensional sums of 
quantum factorials, that are associated to knotted 3-dimensional objects. 
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1.1. The Euler-MacLaurin summation formula. The Euler-MacLaurin summation formula relates 
summation to integration in the following way (see [01 Sec.8]): 

(1) E / (|) = N [ f^ ds + \iSO) /(o)) + R(f, n) 

where the remainder i?(/, A^) has an asymptotic expansion 

(2) R(f,N)~R(f,N) 
in the sense of Poincare, where 

oo 

(3) R(f,N) = E 7^! (/ (2n ~ X) (l) - / (2 ^ 1} (0)) e C[[N-% 

n—l ^ ' 

and B m are the Bernoulli numbers defined by the generating series 

( 4 ) ir^T = E V- 

n— 

Typically, the formal power series R(f,N) is divergent and Gevrey-1, due to the fact that the n-th 
derivative in is not divided by an nl. In the present paper, we discuss an exact form of the Euler- 
MacLaurin summation formula, under a resurgence hypothesis of the function f(x); see Proposition 11.21 
below. 

1.2. Two applications of our exact Euler-MacLaurin summation formula. Our exact form of the 
Euler-MacLaurin summation formula has two applications: in Quantum Topology (where one sometimes 
needs to apply the Euler-MacLaurin summation formula to a resurgent function that has singularities in the 
vertical strip which is perpendicular to the range of summation) , and in Borel summability (with respect to 
e) of difference equations with a small e-parameter. Let us discuss these applications. 

Consider a triple t = (a, b, e) where a, b S Z, b > 0, e = ±1 and the expression 



(5) / t cz) = E« a ^(^ en 

71=0 

where (q) n is the quantum factorial defined by: 

n 

(6) (?)»= 11(1 -<?*), (g) = l. 

fc=i 

Although It(q) does not makes sense when q is inside or outside the unit disk, it does makes sense when 

(a) q is a complex root of unity; in that case It(q) G C. 

(b) q = e 1 /' x ; in that case I t (q) G Q[[l/x]]. 
Given t as above, consider the power series: 

oo 

(7) Lf p (p) = 1 + J2h(e^)p n 

n=0 

(8) L\(p) = B(I t (e^)) 

where B is the Borel transform defined below in Definition 11.11 Our present results, together with |CG3j 
and additional arguments, imply the following theorem, which will be presented in detail in a forthcoming 
publication [CG2] , 

Theorem 1. [CG2] For all i as above, the power series Lf p (p) and if (p) are resurgent functions. 
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In particular, it follows that the generating series Lf p (p) of the Kashaev invariants of the two simplest 
knots, the 3i and 4i (corresponding to t = (0, 1, 1) and t = (2, —1, —1)) are resurgent functions. 

Another application of our exact Euler-MacLaurin formula is to prove parametric resurgence (i.e., resur- 
gence with respect to e) of a formal (WKB) solution to a linear difference equation with a small parameter: 

(9) y(x + e,e) = a(x, e)y(x, e) 
The formal solution of ^ is of the form: 

(10) y(z,e)=exp(if> fc (x)e fc ) . 

\ e k=0 J 

Under suitable hypothesis on a(x, e), Theorem [B] below proves resurgence of the above series for x fixed, and 
constructs an actual solution to © which is asymptotic to the formal solution (flu]) . 

1.3. Known forms of the Euler-MacLaurin summation formula. Before we state our results, let us 
recall what is already known. Suppose that / satisfies the following assumption: 

/ is analytic and satisfies the following bound: 

(11) f(x) = OQx\-) 

for some < 5 < 1 and s > 0, uniformly in the right-half plane 5R(x) > <5. 

For such functions /, Hardy proved in [Hal Sec. 13. 15] that R(f,N) is Borel summable, and that the Borel 
sum agrees with the original sum. In other words, the Borel tranform of R(f, N) can be extended to the ray 
[0, oc), it is integrable of at most exponential growth, and replacing R(f,N) with the corresponding Borel 
sum replaces the asymptotic relation (j2J with an exact identity. 
In a different direction, suppose that 

/ is continuous in the vertical strip < 3J(x) < 1, holomorphic in its interior, and f(x) — o(e 2n ^ x ^) 
as |3(x)| — > oo in the strip, uniformly with respect to 5£(x). 

Then, the Abel-Plana formula states that (see [UJ Sec. 8. 3]): 

( 12 ) £ / (!) = N f Q f(u)du + \{f{l) - /(0)) + i j~ ^^-^ 1 + ^ 2 ;jj^+^ 1 -^ dy . 

1.4. What is a resurgent function? The notion of a resurgent function was introduced and studied by 
Ecalle; see [Eclj . For our purposes, a resurgent function is one that admits endless analytic continuation 
(expect at a countable set of non-accumulating singular points) in the complex plane, and is exponentially 
bounded, that is, satisfies an estimate: 

(13) \f(z)\<Ce a ^ 

for large z. Examples of resurgent functions are meromorphic functions, algebraic functions, or Borel trans- 
forms of solutions of generic differential equations with analytic coefficients. The nth coefficient of the Taylor 
series of a resurgent function around a regular point has a manifest asymptotic expansion with respect to 
1/n that include small exponential corrections; see for example [CG11 Sec. 7]. This property of resurgent 
functions is key in applications to quantum topology, where a main problem is to show the existence of 
asymptotic expansions. For example, an asymptotic expansion of the coefficients of the power series <j3j) is 
almost trivial (for a fixed function /). On the other hand, the existence of asymptotic expansion for the 
coefficients of (x) and (x) (or more generally, i r a ,6,c(e 1/ ' a: )) is a highly non-trivial fact that follows from 
the resurgence of the Borel transform of F a ^ tC {e x / x ); see |CG2| . 

For an introduction to resurgent functions and their properties, we refer the reader to the survey articles 
CNP1, CNP2] [Dl IDP[ IMai ISaj and for a thorough study, the reader may consult Ecalle's original work 
Eel. Ec2 . Let us point out, however, that our main results (Theorems [2] and [4] below) do not require any 
substantial knowledge of resurgence. 
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1.5. Statement of the results. Let us recall a useful definition. 

Definition 1.1. The (formal) Borel transform of a formal power series in 1/N is a formal power series in p 
defined by: 



/ oo \ oo 

(14) B : C^T 1 ]] C[\p]], B E««^TT = E^ 



n=0 / n=0 



-P 



Let Gf{p) denote the Borel transform of the power series R(f, N). 
Theorem 2. If fix) is resurgent and f'{x) is continuous at x = 0, 1 then Gf(p) is given by: 

<"» G > W - i £ ? (/* (' + 5Sl) + /' (' " 5Sl) " f (iHs) - f (-5SS)) 



In particular, G / (p) is resurgent with singularities given by 

(16) TV = {2irinuj, 2irin(LJ — 1) n £ IT ,u — singularity of /'}. 
Let us consider a function /(x) that satisfies the following: 

(Al) / is resurgent with no singularities in the vertical strip < 5R(x) < 1, and f(u) — o(e 27r ' s ' tl ^') as 
|3(it)| — > oo in the strip, uniformly with respect to 

Then, we have the following exact form of the Euler-MacLaurin summation formula. 
Proposition 1.2. Under the hypothesis (Al), for every N 6 N we have: 

N / h \ /"l 1 C°° 

(17) E/ ^ =iV / /(s)ds+ 2 (/(1)_/( ° ))+ / e ~ N " G f(P) d P- 

k— 1 "" 

In particular, the left-hand side of the above equation is the evaluation at N of an analytic function in the 
right hand plane. 

Our proof of Proposition 11.21 allows to generalize to the case that / is resurgent with singularities in the 
vertical strip < 9?(x) < 1; see Theorem|4]in Section 13721 In that case, every singularity A of / in the vertical 
strip gives rise to exponentially small corrections, and the right hand side of Equation (|17p is replaced by a 
transseries. 

Finally, let us give an integral formula for G/(p) which is useful in studying the behavior of Gj(p) for 
large p. 

Theorem 3. With the assumptions of Theorem^ we have: 

<»> <*» - w? r l f^<^ + «--*->) ** 

where 70 is a small circle around oriented counterclockwise. 

Let us end the introduction with some remarks. 

Remark 1.3. Theorem[2j and especially Theorem |4] below provide a new construction of resurgent functions. 
Best known resurgent functions are those that satisfy a difference or differential equation, linear or not; see 
for example [B71 iBrKl ICl] and [Ec2] , 

On the other hand, due to the position and shape of their singularities, the resurgent functions Gf{p) of 
Theorem [2] do not seem to satisfy any differential equations with polynomial coefficients. 

For example, consider the function f(x) = (x — u)~ m where w $ [0, 1] which satisfies the linear differential 
equation with polynomial coefficients: 

(x - w)/'(x) - mf(x) = 0. 

/ is resurgent, with only one singularity at x = w. The corresponding resurgent function Gf(p) of Theorem[5] 
Gf(p) has infinitely many singularities on the rays 27rio;R + , 2iriuMr , 27ri(tj — 1)K + , 2-Ki(to — 1)R~. It seems 
unlikely that Gf(p) satisfies a linear (or a nonlinear) differential equation with polynomial coefficients. 
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Remark 1.4. Let us point out that Theorem [5] implies that the shape of the singularities of Gf(p) is the 
same as that of f'(x). For example, if f'(x) is simply ramified, then so is Gf(p). We recall that a resurgent 
function h(x) is simply ramified if, locally, at each singularity u> of h(x) we have: 

(19) h(x) = P ( — - — | + log(x - u)r(x -ui) + s(x - w) 

\x — lo J 2ni 

where s,r are convergent germs, and P is a polynomial. 

1.6. Acknowledgement. An early version of this paper was presented by the second author at talks in 
Columbia University, Universite Paris VII and Orsay in the spring and fall of 2006. The authors wish 
to thank J. Ecalle for encouraging conversations. The second author wishes to thank G. Masbaum, W. 
Neumann, D. Thurston for their hospitality. 



2. Proof of Theorem [2] 
2.1. Computation of the Borel transform Gf(p). Let © denote the Hadamard product of power series: 

(oo \ / oo \ oo 

n=0 / \n=0 / n=0 

It is classical, and easy to check, that the Hadamard product A®B of two functions A(p) and B(p) analytic 
at p = is also given by an integral formula: 



(21) (A®i?)(p) = -L [ A(s)B(?)d8, 

where 7 is a suitable contour around the origin. For a detailed explanation of the above formula, see [Jul 
p.302] and also [Bol p.245]. 

Let Gf(p) denote the formal Borel transform of the power series in @. Since B m = for odd to > 1, we 
have: 



where 



G,(p) = B 



(l^^-'tD-^-'Ho))^) 

£>2n ( r(2n-l)M\ f(2n-l)tn\\ P 



2n ( f{ 2 n-l)(-i\ _ f (2ri-l)/r,x 

53^(/(™-i)(i)-/(»-i)(o)) *— 

' TO V / TO — 2 

m=2 y 

(£^H«(f:(/""-"w-/'"- i '(o) ; , 

\ T71 — 2 / \ TfL — 2 / 

gi (p) ®92(p) 



(22) fli(p) = E^T 



m! p V eP — 1 n 2 

m=2 * v ^ 
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and 

„m-2 



(23) g 2 (p) = (/^(l)-/^- 13 ^) 



P 



m=2 



(m-2)! 



, / oo 

= i(E(/ < "'- ,, (i)-/ l "'- 1, <o))(£_ 1)!/ 

= ^ W + P) - /(P) -/(D + /(°» = /(I +P) " /»■ 

Consider positive numbers r and 8 such that <?i(p) is analytic for \p\ < r (eg, r < 27r) and 32 (p) is analytic 
for |p| < d-the latter is possible by Equation (f2"5| and our assumptions on /. 
Now, Equation (f2"Tj) implies that 

1 /" . . /p\ 



,m — 1 



(24) G ^ = ^iJ^Ml) 

for all p with \p\ < Sr, where 7 is a circle around with radius r < r . 

With our assumptions, when \p\ < Sr, the function s 1— > g2{p/s) has no singularities outside of 7. Thus, 
outside of 7, the singularities of gi(s)g2(p/ s)/s are simple poles at the points 2mn for n € Z*, with residues 

R es ( 5i( s )ff2 (-) -, s = 27rin ) = Res(g 1 (s), s = 2mn) g 2 ( 77^) tt^" 
\ \s/ s J \lmn) Zmn 

1 2 f P 1 

Moreover, <?i(s) = 0(1/ s) when the distance of s from 2-kiL* is greater than 0.1 and gi(p/s) = 0(1) for s 
large, thus the integrand vanishes at infinity. 

We now enlarge the circle 7 and collect the corresponding residues by Cauchy's theorem. Using the above 
calculation of the residue and Equation (|23j) . it follows that 



1 00 1 

n=l 

= lyl f 1 + _£J) + /' f 1 - - /' (JL-} r (~^- 

47r 2 n 2 V V 27rm/ V 2wins \2mnJ V 2?rin 
n=i 

Since /'(x) is regular at x = 0, 1, it follows that the above series is convergent for p G C — TV, where TV is 
defined in (|16[) . In addition, we conclude that Gf(p) has endless analytic continuation with singularities in 
TV. 

It remains to prove that Gf(p) is exponentially bounded, assuming that / is. If / is exponentially bounded, 
Cauchy's formula implies that /' is exponentially bounded. Then, we have: 



''( 1+ 24)l^ Cexp ( a 



2nin 



Cr 1 cxp I 



Thus, 



CV + 1) C(e" + 1) / |p|\ 

|G ' WI * 6XP { a 2^ ) g = 6XP v a 27 J ' 

This completes the proof of Theorem [2j □ 
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3. An exact form of Euler-Maclaurin summation formula 

3.1. Proof of Proposition [T72l Proposition 11.21 follows easily from the Abel-Plana formula; see Appendix 
IA1 However, we give a proof of Proposition 1 1 . 2 1 that allows us to generalize to Theorem [4] below. 

Consider a resurgent function / that satisfies the assumptions (Al), and let us introduce the function 

h(u) = —f(u) 777 rvp- 

and the contour Tr.$ which is a rectangle oriented counterclockwise with vertices — iR, 1 — iR, 1 + iR, iR 
that excludes the points 0, 1 together with small semicircles of radius 6 at the points and 1. 



1 



C, 



Figure 1. The contours Ci , Ci , C3 , C4 of the critical strip, and a truncated contour Fr,s 

Due to our assumptions on /, the singularities of h(u) inside Tr.$ are simple poles at k/N with residue 
f(k/N) I (27ri) for k = 1, . . . , N — 1. The residue theorem implies that 



(25) 



k=l 



du. 



Let s (resp. T R s ) denote the upper (resp. lower) part of the contour T. Since f(x) has no singularities 
in 5R(u) G [0, 1], the residue theorem implies that 



(26) 



N / f(u)du 



N 



10 2 J r + 

Adding up Equations (f2"5")) , (|2"5)) and using 

1 z + r 1 



N f 

f{u)du - — I f(u)du. 



2z- z- 1 
2Z-Z- 1 



1 



l-z~ a 
1 

z 2 - 1 



we obtain that 



(27) £ / (£) - * Si = * / r , , r^L^ + ^4 , ^n d " 

Now let i? — > 00. Due to assumption (Al), the integrals over the horizontal parts of s approach zero. 
Next, let S — > 0. Since / is continuous, the integral around the quarter circle that links S to id tends to 
— /(0)/4. The other quarter circles are treated similarly. 

Thus, we have: 
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f(u)du=-(f(l)-f(0)) 



/(«)-/(!) 



C-2 



1-e 



-2-rriNu 



du 



m - /(o) 



c 4 



1 -e 



-2-KiNu 



du 



/(«)-/(!) 



27riNu 



du 



/(«) - /(o) 



c 3 



du. 



Consider now the corresponding function Gfijp) from Theorem [5J We have: 

G/(p) = Gi(p) + G a (p) + Ga(p) + G 4 (p) 

where 

oo 



Gi(p) 
G a (p) 
G 3 (p) 
G 4 (p) 



— T-f fi 

4tt 2 ^ n 2 J \ 

n—1 
1 00 i 

— T-f fi 

4^2 n 2J V 

71=1 

-— F-z'f— 

n=l 
1 °° 1 
^2 T1 2 J 



P 



47r 2 — ' n 

71 = 1 



Inn 
P 



Consider the contours C\, C2, C3, C4 on the boundary of our strip, as shown in Figure [TJ 
We claim that the Laplace transform of the Gj (p) for j = 1, . . . , 4 is given by: 



(28) 



- Np G j (p) = < 



/(«)-/(!) 



c 5 - 





- 1 


/(«)- 


/(O) 






fin)- 


/(I) 


^ _ p — 2-rtiNu 


/(«)" 


/(0) 



for j = 1 
du for j = 3 
e?u for j = 2 
, ., du for 7 = 4. 

2mNu J 



Let us show this for j = 3; the other integrals are treated in the same way. We compute as follows: 



/>oo 1 />oo 00 -1 



P 



2tH /- E 



00 e -2niNnu 



2nin 
f(u)du 



dp 



C 3 . 



-2ixiNu\ j?/ 



)f{u)du 



= N 



•c 3 

f{u) - /(0) 

a 2niNu ^ 



du 



by interchanging sum and integral 
hyp = 2mnu 
by ([29|) 

by integration by parts 



where 
(29) 



^ e -2niNnu 



E 



log(l - e 



-2-KiNu 



This concludes the proof of Proposition 1 1.21 in case / satisfies (Al). 



□ 
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Let us end this section with a remark. 

Remark 3.1. If f(x) = e cx , one may verify Equation (|15p directly by using the Mittag-Lefner decomposition 
of the function x/(e x — 1). 

3.2. Euler-MacLaurin summation for functions with singularities in the vertical strip. In this 
section we consider a function / that satisfies the following assumptions: 

(A2) / is resurgent, and let A denote its set of singularities on the critical strip < $t(%) < 1. We assume 
that A £ [0,1] for all A G A, and f(u) = o(e 27r l 3 ^>l) as |3r(u)| — ► co in the strip, uniformly with 
respect to SR(u). We also assume that on a vertical ray A + iK + , we have /(u)e _27r ^^ u - ) G L 1 (A + zE + ). 

(A3) For every A € A there exist a holomorphic germ h\(u) and real numbers a\, (3\ so that for u near 
we have: 

(30) f {u + X) = u a >- (log uf^h x (u). 

(A4) For simplicity, let us also assume that 5ft(A) ^ 5R(A') for A ^ A', and that A is a finite set. 
Let 

(31) (£G)(x) = / e~ xp G(p)dp 



denote the Laplace transform of G(p). We denote by f\(u) the variation (or jump) of the multivalued 
function /(it + A) at it; where u lies the vertical ray starting at (see for example, [Ma]). We also define: 

(32) G fXm (p)=i-^—h' "' 



27rm \27rm 

In case f(u + A) is single- valued then Gf t \ t m(p) is a distribution supported at p = 0. 
Then, we have the following exact form of the Euler-MacLaurin summation formula. 

Theorem 4. (a) If / satisfies (A1-A4) and a\ > — 1 for all A G A, then for every f GNwe have: 
(33) = iV^ 1 /(^ + 5(/(l)-/(0)) + (£G / )(iV) 



DC 



+ Ne 2mXN E e 2 ™ AmW (£G /Am )(A0 

X:9f(A)>0m=0 
oo 

+ Ne -^XN ^ ^ e- 27rlAmAr (£G / , A , m )(A^) 

A:9(A)<0"i=0 

(b) If some ax < — 1, Equation (|33|) is true after integration by parts M-times where M > niaxj, :ai <-i[- oja]. 

Proof. Without loss of generality, let us assume that / has a single singularity A in the vertical strip < 
3?(x) < 1 with 3(A) > 0. 

Use the modified contour Trs \ in Figured] 

Let Hr denote the portion of V r g x that consists of the truncated Hankel contour around A, and Sr j( 5 = 
r,R,c5,A — Hr. Equations ([23]) and (|2"b]l become: 

AT r niNu I -iriNu " / U \ AT C e -niNu i —TtiNu 

(34) -v/ /(") w w ^ + =v / /(") w w <*" 



Hr k=l 



and 



AT /■ f 1 TV /" AT /" 

- / /( U )d u - A / f{u)du = - / /(«)d« - - / 
1 Jh r Jo z Js+z 1 Js\ 



(35) - — / f{u)du - N / f(u)du = — I f{u)du - — I f{u)du. 

Adding up, the extra contribution from Hr becomes: 
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iR 




(?) 



Figure 2. The modified truncated contour Tr^^x on the left and a Hankel contour H on the right. 



(36) ~n[ f } U l lNu du 

Jh r 1 - e 

Now let R — > oo. Notice that f(u + A) is uniformly L 1 for u near iff a\ > —1 for all A. Using this and our 
integrability assumption (A2), it follows that in the limit the above integral equals to 

IX = - N j H l -{-2«iNu dU 

= -^h — f f s L + ^ s . 

Now, A = Ai + i\2 with A2 > 0, and we may write 

1 1 



^ _ e -27iiN(\+is) I _ fjj—N e 2irNs 



00 



N -2-rrNs Nm -2Nms 



27rJVs 

m=0 

where o> = e 27 ™ A satishcs M < 1. Thus, 



-iVu;* ]T w «™(£G A , m )(JV). 



m=0 



Part (a) of Theorem |4] follows. Part (b) follows from the fact that if f(u + A) has a local expansion of the 
form (|30|) . and F^(u) = f(u), then f(u + A) as a local expansion of the form: 

(37) F(u + A) = u a ^ +s (\oguY^H x (u) 

for a holomorphic germ H\(u). Cf. also CI, Thm.l]. □ 

3.3. Euler-MacLaurin with logarithmic singularities at x = 0. In this section we consider functions 
f(x) that have a logarithmic singularity at x — 0. Motivated by our applications to quantum topology, we 
consider functions / of the form: 



(38) f(x) = clogx + g(x) 

where g that satisfies (Al), and c £ C. Let us define 
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(39) H(p) = -{-^- r -l + ^ 



p 2 \eP-l 2 / 

It is easy to see that H(p) is analytic at p = 0. In fact, the Taylor series of H at p — is given by: 



(40) S^ 2 "" 
Theorem 5. Under the above hypothesis, for every N £ N we have: 

( 41 ) T,Hm)= N f(s)ds+-logN+-log(2n) + -(g(l)-g(0)) + C(G g +cH)(N). 

fc=l V / "/0 

Proof. Since / is given by (|38p . we have: 



V A / ' B \N N J ^ a \N 

k=l v 7 v 7 k=l v 

Recall now from |Ha|, Sec. 13. 15] the following exact form of Stirling's formula: 

( iV \ i l 

(42) l0g (iV^J =2 l0 S Ar - iV +2 l0g(27r) + ( ' CjH " )(iV) - 

Applying Proposition 1 1 . 21 to g gives: 

JV 



f k \ I" 1 1 

E^ lv =iV / 9(^+2(-9( 1 )-5(0)) + (^G 9 )(iV). 

Adding up, and using 

N f f(s)ds = N [ g(s)ds + Nc[ log sds = N [ g(s)ds - Nc 
Jo Jo Jo Jo 

we obtain (|4T1). The result follows. □ 



4. Parametric resurgence of difference equations with a parameter 
Consider the first order linear difference equation with a small parameter e: 

(43) y(x + e, e) = a(x, e)y(x, e) 

where a(x, e) is smooth. (|43|) has a unique formal solution (often called a WKB solution) of the form: 

(44) y(x,e) = e-^o 

where £}(0) = 0. See for example, |CCj and |GGj . For simplicity, suppose that a(x, e) = a(x) is independent 
of e. Under the stated assumptions, the next theorem gives an exact solution to (|43|) which is asymptotic to 
the formal solution (|4"4")l . 

Theorem 6. (a) For all x such that s — ► log a(sx) satisfies (Al) we have: 

(45) _ V^ fc (x)e fc - - / log a(q)dq-- log a(x) + -logo(0) + / e-«/ e G( g , a;)dg 
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where 



(46) G ( q , A e 4 f fl ? + t} + a '\ x ~ tI - 4#4 - 4^441 

n—l \ V 2irtn ) \ 27Ttn / \2-KinJ V 2-jnn / / 

1 f°° f u ( log a(s). pu m N log a(x + s), eh eh , . 



(27ri) 3 y J 7o e u - 1 V s 2 (x + s) 

where 70 is a small circle around oriented counterclockwise 
(b) Moreover, (|4"5|) has a solution y(x, e) of the form: 



(47) yfa e ) = ^2M exp (1 £ loga(q)dq + J™ e~^G(q, x)dq^j . 

Remark 4.1. It follows that the singularities of G(q, x) are of the form 2itin\ or 2irin( A — x) where b£Z* and 
A is a singularity of log a. These type of singularities appear in parametric (i.e., co-equational) resurgence 
of Ecalle; see |Ec2j . 

The proof of Theorem [6] indicates the close relation between the Euler-MacLaurin summation formula 
and the formal solutions of a linear difference equation with a parameter. 

From that point of view, resurgence of Gf(p) translates to parametric resurgence of formal solutions of 
linear difference equations. In the case of formal solutions of linear differential equations with a parameter, 
Ecalle shows that their singularities are of the form n{a.i — x) for n = —1, 1, 2, 3, ... ; see |Ec2[ Eqn.(6.9)]. 

Proof, (a) Let z(x, e) = logy(x,e). Taking the logarithm of (j43|) . it follows that 

z(ke + e, e) = loga(fce) + z(ke, e). 
Summing up for k = 0, . . . , N — 1 and using the variable 

(48) x = Ne 
we obtain that: 



7V-1 



z(x, e) - z(0, e) = ^ loga(fce) 



fc=0 

N 



loga(x ) +loga(0) +^2\og a(xk/N). 



fc=i 



Let us fix x and apply Proposition 1 1 . 21 to the function s — > log a(sx). We obtain that 



z(x,e) - z(0,e) = N log a(xs)ds -- log a(x) + - log a(0) + / e~ Np H(p,x)dp 



1 f x 1 1 f°° 

- I log a(s)ds - -log a(x) + -log a(0) + / e~ xp/e H(p,x)dp 

6 Jo 2 2 Jo 



e 

where by Theorems [2] and [3] we have: 



f log a(s)ds-- log a(x) + - log a(0)+ / e^'H (^-,x) 
2 2 J \x I 



1 rr « ^iog«H (e ^ +e _ lfe) _ logaqi +,)x) (e ^_ +e _^_A ^ 



(2ttz) 3 7 7 7n e" - 1 V s 2 (1 + s) 2 
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Since H(q/x, x)/x — G(q, x) (where G(q, x) is given by (|3S]0 and z(0, e) = 0, it follows that for all e > and 
k G Z we have: 

(49) z(ke + e, e) = loga(fce) + z(ke, e). 

To prove (b), let us consider the difference 

E(x,e) = e(logy(x + e, e) - loga(ir) - \ogy(x,e)). 
It follows by definition that E(x,e) is analytic in (x, e). Thus, 

oo 

J5(a:,e) = cyjcV. 
*,j=o 

Moreover, (|4"9"| implies that for all e > and all k € Z we have: 

oo 

= F(fce,e) = 2 r, ; /,•','••'. 

Thus Cjj = for all i, j and E(x, e) = 0. This completes the proof of (b). 

The definition of y(x, e) by a Laplace integral and Watson's lemma (see [Oj Sec. 4. 3.1]) implies that 



l(x,e) ~ -^T<p k (x)e k 



e 

k=0 



for analytic functions 4> k (x) that satisfy <pk (0) = 0. Since a formal WKB solution given by (|44|) is unique, it 
follows that Ff.(x) — 4>k{x) for all fc. Thus, (a) follows. □ 

Remark 4.2. Theorem [5] can be generalized when 

oo 

a(x,e) = ^a k (x)e k 

k=0 

is analytic with respect to (x, e), and the coefficients a,k{x) are resurgent functions. It may also be generalized 
to the case of higher order linear difference equations with a parameter. This will be explained elsewhere. 

Remark 4.3. The reader may compare Theorem [S] with the results of the last section of |Sa| . 

5. AN INTEGRAL FORMULA FOR Gf(p) 

In this section we give a proof of Theorem [3l We follow the ideas of |C2| to convert the sum of Equation 



(fT5]l into an integral. Let us show that 
(50) 



£-[n 2j \2ninJ 2ni J J l0 s 2 (e u - 1) 



and similarly for the sum of the other three terms in (fT"5)) . 

To prove Equation (|50p . we first expand /' at p — 0, then take a Laplace transform with respect to the 
summation variable n, interchange the order of summation and sum the geometric series. We obtain that: 



v^i.m yy / (j+1) (0) /PV l 
^ n 2 \2mn) j\ \2iri) n'J+ 2 



71=1 j=0 



oo oo 



fjf^U ■ l)!(2*W 

f^Jo j!(j + l)!(2«y e*-l 



14 



OVIDIU COSTIN AND STAVROS GAROUFALIDIS 



Using Cauchy's formula 

/O'+i)(0) _ 1 f f( s ) 



- ds 
(7 + 1)! 2 m J 10 s^ ' 



and interchanging summation and integration it follows that 

f(s)p> u j+1 1 

3=0 



f 1 P \ = 1 y f f /W 1 
^ n 2 V27rin/ 2tti ^ 7 7 j!sJ+ 2 (27ri>? e u - 1 



^70 

i_ /-°° /■ uf{ s ) Y-UP u y dcdu 

2ni J J l0 s 2 (e« - 1) £j j! V2tuJ 

i r /■ «/w 



v-i / " _ "n 

70 



Wo y 70 s 2 (e«-l) 



The interchanges of summation and integration are justified by dominated convergence. This concludes the 
proof of ([50)) and Theorem [3] □ 

Appendix A. 

For completeness, let us show how the Abel-Plana formula implies Proposition 1 1.21 With the notation as 
in Proposition II. 21 we claim that for every N € N we have: 



Adding up, and using the Abel-Plana formula (jl2|l . gives a proof of Proposition 11.21 Let us give the proof 
of (f5Tj) and leave the rest as an exercise. For y > 0, we have e~ 2vy < 1 and the geometric series gives: 



(55) -= = y e - 2 ™^. 

n— 1 

Interchanging summation and integration, changing variables 27rny = Np and integrating by parts (justified 
by the hypothesis (Al)), we obtain that 



n— 1 u 

- 5 Eire-*)' (/('-55s) -/<•>)* 

n— 1 u 

4tt 2 ^ n 2 Jr. J \ 2mnJ 1 

n— 1 u 
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